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Abstract 

A recently developed dynamical mean-field theory in the iterated per- 
turbation theory approximation was used as a basis for construction of 
the "first principles" calculation scheme for investigating electronic struc- 
ture of strongly correlated electron systems. This scheme is based on 
Local Density Approximation (LDA) in the framework of the Linearized 
Muffin- Tin-Orbitals (LMTO) method. The classical example of the doped 
Mott-insulator Lai-^Sr^TiOa was studied by the new method and the 
results showed qualitative improvement in agreement with experimental 
photoemission spectra. 

1 Introduction 

The accurate calculation of the electronic structure of materials starting from 
first principles is a challenging problem in condensed matter science since un- 
fortunately, except for small molecules, it is impossible to solve many-electron 
problem without severe approximations. 

For materials where the kinetic energy of the electrons is more important 
than the Coulomb interactions, the most successful first principles method is 
the Density Functional theory (DFT) within the Local (Spin-) Density Ap- 
proximation (L(S)DA)Q, where the many-body problem is mapped into a non- 
interacting system with a one-electron exchange-correlation potential approxi- 
mated by that of the homogeneous electron gas. 

It is by now, generally accepted that the spin density functional theory in 
the local approximation is a reliable starting point for first principle calculations 
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of material properties of weakly correlated solids (For a review see pj). The 
situation is very different when we consider more strongly correlated materials, 
(systems containing f and d electrons). In a very simplified view LDA can be 
regarded as a Hartree-Fock approximation with orbital-independent (averaged) 
one-electron potential. This approximation is very crude for strongly correlated 
systems, where the on-cite Coulomb interaction between d- (or f-) electrons of 
transition metal (or rare-earth metal) ions (Coulomb parameter U) is strong 
enough to overcome kinetic energy which is of the order of band width W . In 
the result LDA gives qualitatively wrong answer even for such simple systems 
as Mott insulators with integer number of electrons per cite (so-called " undoped 
Mott insulators") . For example insulators CoO and La2Cu04 are predicted to 
be metallic by LDA. 

The search for a "first principle" computational scheme of physical proper- 
ties of strongly correlated electron systems which is as successful as the LDA 
in weakly correlated systems, is highly desirable given the considerable impor- 
tance of this class of materials and is a subject of intensive current research. 
Notable examples of first principle schemes that have been applied to srongly 
correlated electron systems are the LDA+U method S which is akin to orbital- 
spin-unrestricted Hartree-Fock method using a basis of LDA wave functions, ab 
initio unrestricted Hartree Fock calculations B and the use of constrained LDA 
to derive model parameters of model hamiltonians which are then treated by 
exact diagonalization of small clusters or other approximations . 

Many interesting effects, such as orbital and charge ordering in transition 
metal compounds were successfully described by LDA+U method || . However 
for strongly correlated metals Hartree-Fock approximation is too crude and more 
sophisticated approaches are needed. 

Recently the dynamical mean-field theory was developed B which is based 
on the mapping of lattice models onto quantum impurity models subject to a 
self-consistency condition. The resulting impurity model can be solved by var- 
ious approaches (Quantum Monte Carlo, exact diagonalization) but the most 
promising for the possible use in "realistic" calculation scheme is Iterated Per- 
turbation Theory (IPT) approximation, which was proved to give results in a 
good agreement with more rigorous methods. 

This paper is the first in a series where we plan to integrate recent devel- 
ompements of the dynamical mean field approach with state of the art band 
structure calculation techniques to generate an "ab initio" scheme for the cal- 
culation of the electronic structure of correlated solids. For a review of the 
historical development of the dynamical mean field approach in its various im- 
plementations see ref . In this paper we implement the dynamical mean-field 
theory in the iterated perturbation theory approximation, and carry out the 
band structure calculations using a LMTO basis . The calculational scheme is 
described in section 2 . We present results obtained applying this method to 
Lai_ ;E Sr x Ti03 which is a classical example of strongly correlated metal. 
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2 The calculation scheme 



In order to be able to implement the achievements of Hubbard model 
theory to LDA one needs the method which could be mapped on tight-binding 
model. The Linearized Muffin-Tin Orbitals (LMTO) method in orthogonal 
approximation § can be naturally presented as tight-binding calculation scheme 
(in real space representation): 

HlmTO — i^ilmjl'm' e iina m (T +tilm,jl'rn'%l m a^jl' 1 rn'(r) (1) 

(i - site index, Im - orbital indexes). 

As we have mentioned above, LDA one-electron potential is orbital - inde- 
pendent and Coulomb interaction between d-electrons is taken into account in 
this potential in an averaged way. In order to generalize this Hamiltonian by 
including Coulomb correlations, one must add interaction term: 

H-int 2 ^ ^ Tlii 7na Tln m f a ' (2) 

ilrara' gg' mG^ra' g' 

We neglected for a while exchange terms and dependence of Coulomb pa- 
rameter U on the particular pair of orbitals mm' . Through the following we will 
assume that only for one shell Id of one type of atoms id (for example d-orbitals 
of the transition metal ions) Coulomb interaction will be taken into account 
(Uu — USu t i d i d ), and in the following iesndex il will be omitted. All other or- 
bitals will be considered as resulting in the itinerant bands and well described 
by LDA. Such separation of the electronic states into localized and itinerant is 
close in spirit to the Anderson model. 

To avoid double-counting one must in the same time subtract the averaged 
Coulomb interaction energy term, which we assume is present in LDA. Un- 
fortunately there is no direct connection between Hubbard model and LDA 
(because LDA is based on the homogeneous electron gas theory and not on the 
localized atomic type orbitals representation) and it is impossible to express 
rigorously LDA-energy through d-d Coulomb interaction parameter U. How- 
ever it is known that LDA total energy as a function of the total number of 
electrons is a good approximation and the value of the Coulomb parameter U 
obtained in LDA calculation agrees well with experimental data and results of 
the more rigorous calculations ||.That leads us to the suggestion that a good 
approximation for the LDA part of the Coulomb interaction energy will be: 

Ecoui = 7;Un d (n d - 1) (3) 
( n d = X) n ma total number of d-electrons). 
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In LDA-Hamiltonian e d has a meaning of the LDA-one-electron eigenvalue 
for d-orbitals. It is known that in LDA eigenvalue is the derivative of the total 
energy over the occupancy of the orbital: 

e d = 4~ Elda ( 4 ) 
dn d 

If we want to introduce new e d where d-d Coulomb interaction is excluded 
we must define them as: 

e d =^- d {E LDA -Ecoui) = e d -U{n d - 1 -) (5) 

Then new Hamiltonian will have the form: 

H = H + Hint 

H ^ ^ (^ilnijl' m' ^il n ilmcr tilnijl' ni' Cii ma Cjl' rn'a^) (6) 

ilmjjl'm' .a 

In reciprocal space matrix elements of the operator H° are: 

Hqlm,q'l>m'0 S ) = ^qh^q'Vra'^) ~ ^qlm,q'l'm' 5ql,i d l d U(n d - -) (7) 

(q is an index of the atom in the elementary unit cell). 

In the dynamical mean-field theory the effect of Coulomb correlation is de- 
scribed by self-energy operator in local approximation. The Green function 
is: 

G q lm,q>l>,n>{iu) = ^- J rfk [w + [i - H° lm q , l>m ,(k) - Sqlm^l'm'SqLi^iiu))}- 1 

(8) 

([...] 1 means inversion of the matrix, integration is over Brillouin zone, /i is 
chemical potential, Vb is a volume of Brillouin zone). 

In the following we will consider paramagnetic case, orbital and spin degen- 
erate system, so that self-energy S(io;) does not depend on orbital and spin 
indexes. One can define effective Anderson model Green function through: 

G(iu) = G id i dm , id i dm (iuj) = (iuj + fi- A(iu) - E(w)) -1 (9) 

where A(iw) is effective impurity hybridization function. The effective medium 
"bath" Green function G° is defined as: 

G°(iu>) = (iw + Jl- A(iw))" 1 = (G _1 (iw) + S(iw) +J1- /i)" 1 (10) 

(jl is chemical potential of the effective medium). 
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The chemical potential of the effective medium \x is varied to satisfy Luttinger 
theorem condition: 

lj2e^ + G(i. n ) 1 ^E(i. n ) = (11) 

In iterated perturbation theory approximation the anzatz for the self-energy 
is based on the second order perturbation theory term calculated with "bath" 
Green function G°: 



Z°(iw s ) = -(7V-l)f/ 2 -i^^G (z W „ l + z J5n )G ^ m )G ^ s -zp n ) (12) 
N is a degeneracy of orbitals including spin, (3 = Matsubara frequencies 



w s = ( 2s + 1 ) 7r -p Ti = integer numbers. 

The term E° is renormalizcd to insure correct atomic limit: 



1 - BYP(iu) 



E(to) = ETn(jy-l) + - 5^ (13) 



(n is orbital occupation n = 4 ^ e " J ™ 0+ G'(icj n )), 



5 = + g 

t/ 2 (^~i)n (i-n ) 



A = n[l - (jV - M + (N- 2)D[n] 
n (l-n ) 

correlation function D[n] =< nn >cpa is calculated using Coherent Poten- 
tial Approximation (CPA) for the Green function with parameter Sfj, chosen to 
preserve orbital occupation n : 

[l-n(JV-l)] n(JV-l) 

GCPA^U) - + ■ t-t— „ (17) 

+ [I — £\{IUI) +0(1 ILO + fl — /\(IU1) — U + Of! 

lj2e^ 0+ GcPA(ioJ n ) (18) 



p 



tLU-,, 



D[n] = «ye ,u - 0+ : — (19) 

The Matsubara frequency convolution in (|l2l ) was performed with time vari- 
ables representation using Fast Fourier Transform algorithm for transition from 
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energy to time variables and back: 



G V) = ^E e ~"" TG °(^) ( 2 °) 

iuj n 

E(r) = — (iV - 1)U 2 G°(t)G°(t)G°(—t) (21) 

P 

E°(iw n ) = /dre M " r E(T) (22) 



The serious problem is to perform integration in k-space over Brillouin zone. 
For this we used generalized Lambin-Vigneron algorithm [[L0| . We define new 
matrix H(k, z) as: 

ff(k,«) = H°(k)+E(z) (23) 

z - complex energy, term E(z) is added only to diagonal element of //-matrix 
corresponding to d-orbitals. In this matrix notations Green function is : 

G{z) = ±- J dklz-H^z)]- 1 (24) 

After diagonalization, H (k, z) matrix can be expressed through diagonal matrix 
of its eigenvalues D(k, z) and eigenvectors matrix f/(k, z) : 

H(k, z) = U{k, z)D(k, z)U-\k, z) (25) 

and Green function: 

G{z) = ^r [ dkU(k, z)[z - D(k, z^U^ik, z) (26) 
Vb J 

A particular matrix element of Green function is calculated as: 

^ 1 f UiJk,z)U-}(k,z) 

In analytical tetrahedron method the irreducible wedge of the Brillouin zone 
is divided into a set of tetrahedra and the total integral is calculated as a sum 
over the tetrahedra. To perform integration over a given tetrahedron with four 
corners at vectors k^ (i = 1, 2, 3, 4) denominator of the fraction in equation ( p7| ) 
is interpolated as a linear function in k-space. In the result the integral over 
one tetrahedron is expressed through the values of numerator and denominator 
at the corners of the tetrahedron: 

1 f U m (k,z)U^(k,z) * „ 

n \ f / n i — 1 



G 



V is tetrahedron volume 

„ _ (z - D n (ki,z)) 2 



11 {D n {1t k ,z)-D n Qsi,z)) 



(29) 



(j - £> n (kj, z)) 3 ln[(z - £> n (kj, - D n (kj, z)] 

^ II (£>n(k*, «)-£>„ (k Js z)) (I?„(k 4 ,z)-I?„(k,,z)) 

The self-energy £(iu; n ) and Green function G(iu) n ) are calculated at the 
imaginary Matsubara frequencies iu> n — iir(2n+ l)//3 . It is enough to calculate 
expectation values, such as orbital occupancies n , but in order to calculate 
spectral properties one need to know Green function on the real axis. The 
real axis equivalent of equations (|l^) is much more complicated and hard to 
implement numerically than Matsubara frequencies version. It is much more 
convenient to perform analytical continuation from imaginary energy values to 
the real ones. For such continuation we have used Pade approximant algorithm 
[jllj. If one has a set of the complex energies Zi (i — l,...,M) and the set 
of values of the analytical function m , then the approximant is defined as 
continued fraction: 

n / a o,x ai{z ~ z 2) a M (z-z M -i) 

Cm(z) = — — ... [60) 

where the coefficients are to be determined so that: 

C M (z i )=u ii i = l,...,M (31) 
The coefficients aj are then given by the recursion: 

a t = g%(zi), gifa) — Ui, i — 1, M (32) 

f \ 9p-i(z P -i) -g p -i(z) 

g P {z) = —, \ r^— , P>2 33 

(z - z p -i)g p ^i{z) 

The recursion formula for continued fraction finally yields: 

C M (z) = Am(z)/B m (z) (34) 



where 



and 



A n+1 (z) = A n (z) + (z - z n )a n+ iA n -i(z) 

B n +i(z) = B n (z) + {z - z n )a n+1 B n _i(z) (35) 

A = 0, Ax = oi, B = B 1 = 1 



We have found that the most convenient way is to use analytical continuation 
not for the Green function G but only for self-energy E, and then to calculate 
G directly on the real axis through the Brillouin zone integration (12^ 
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3 Results 



We have applied the above described calculation scheme to the doped Mott 
insulator Lai_ x Sr :z .Ti03. LaTi03 is a Pauli-paramagnetic metal at room tem- 
perature and below TV = 125 K antifcrromagnetic insulator with a very small 
gap value ( 0.2 eV) . Doping by a very small value of Sr (few percent) leads to the 
transition to paramagnetic metal with a large effective mass. As photoemission 
spectra of this system also show strong deviation from the noninteracting elec- 
trons picture, Lai-^Sr^TiOa is regarded as an example of strongly correlated 
metal. 

The crystal structure of LaTiC>3 is slightly distorted cubic perovskite. The Ti 
ions have octahedral coordination of oxygen ions and ti g -e g crystal field splitting 
of d-shell is strong enough to survive in solid. On Fig.l the total and partial 
DOS of paramagnetic LaTi03 are presented as obtained in LDA calculations 
(LMTO method). On 3 eV above 02p-band there is Ti-3d-band splitted on i 2 g 
and e g subbands which are well separated from each other. Ti 4+ -ions have d 1 
configuration and ti g band is 1/6 filled. 

As only ti g band is partially filled and e g band is completely empty, it is 
reasonable to consider Coulomb correlations between ti g — electrons only and 



degeneracy factor N in Eq. (12) is equal 6. The value of Coulomb parameter 
U was calculated by the supercell procedure Q regarding only ti g — electrons as 
localized ones and allowing e g —electrons participate in the screening. This cal- 
culation resulted in a value 3 eV. As the localization must lead to the energy gap 
between electrons with the same spin, the effective Coulomb interaction will be 
reduced by the value of exchange parameter J=l eV. So we have used effective 
Coulomb parameter U e ff=2 eV. The results of the calculation for x=0.06 (dop- 
ing by Sr was immitated by the decreasing on x the total number of electrons) 
are presented in the form of the i2g-DOS on Fig. 2. Its general form is the same 
as for model calculations: strong quasiparticle peak on the Fermi energy and 
incoherent subbands below and above it corresponding to the lower and upper 
Hubbard bands. 

The appearance of the incoherent lower Hubbard band in our DOS leads to 
qualitatively better agreement with photoemission spectra. On Fig. 3 the exper- 
imental XPS for Lai^^Sr^TiOs (x=0.06) jl^] is presented with non-interacting 
(LDA) and interacting (IPT) occupied DOS broadened to imitate experimental 
resolution. The main correlation effect: simultaneous presence of coherent and 
incoherent band in XPS is successfully reproduced in IPT calculation. However, 
as one can see, IPT overestimates the strength of the coherent subband. 

4 Conclusions 

In this publication we described how one can interface methods for realistic 
band structure calculations with the recently developed dynamical mean field 
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technique to obtain a fully "ab initio" method for calculating the electronic 
spectra of solids. 

With respect to earlier calculations, this work introduces several method- 
ological advances: the dynamical mean field equations are incorporated into 
a realistic electronic structure calculation scheme, with parameters obtained 
from a first principle calculation and with the realistic orbital degeneracy of the 
compound. 

To check our method we applied to doped titanates for which a large body 
of model calculation studies using dynamical mean field theory exists. There 
results are very encouraging considering the experimental uncertainties of the 
analysis of the photoemission spectra of these compounds. 

We have used two relative accurate (but still approximate ) methods for 
the solution of the band structure aspect and the correlation aspects of this 
problem: the LMTO in the ASA approximation and the IPT approximation. 
In principle, one can use other techniques for handling these two aspects of the 
problem and further application to more complicated materials are necessary to 
determine the degree of quantitative accuracy of the method. 
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5 Figure captions 

Fig. 1. Noninteracting (U=0) total and partial density of states (DOS) for 
LaTi0 3 . 

Fig. 2. Partial fa g ) DOS obtained in IPT calculations in comparison with 
noninteracting DOS. 

Fig. 3. Experimental and theoretical photoemission spectra of Lai-^Sr^TiOa 
(x=0.06). 
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Fig. 2 
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